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ABSTRACT 

Factors of importance in power plant site selection for 
once - through cooling water systems situated on sea coasts 
are reviewed. Principles of recirculation management are in- 
dicated and briefly discussed in view of the implications of 
this approach on the environmental aspects of the problem. 

The need to recognize that the waste heat problem is not sol- 
ved by prescribing zero recirculation is stressed. 


INTRODUCTION 

Selection of coastal sites for large power plants must consi- 
der the recirculation aspects as well as the possibility of 
environmental side effects. 

All heat rejected from power plants is ultimately dissipa- 
ted to the atmosphere which is achieved by raising the wa- 
ter temperatures above their natural ambient level over a 
substantial area. Hence, for a given site and heat load the 
excess temperature distribution can only to a limited ex- 
tent be managed by engineering design. 

Recirculation resulting from the rise of ambient temperatu- 
re in the vicinity of the intake can be made subject to a 
rational management. Investments in marine installations 
should be adjusted, not necessarily to minimize recircula- 
tion, but to minimize the total cost of power loss due to 
recirculation and costs of marine installations. 

An approach based on recirculation management principles 
does not appear to be in conflict with considerations with 
respect to the question of thermal impact on the acquatic 
environment at a given site. 
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MANAGEMENT PRINCIPLES 

In general terms the objective of waste heat management is 
to minimize the sum of construction and plant operation 
costs and possible damages to the equatic invironment. Con- 
sequently a possible objective function is: 

F = R + D + C (1) 

where R is the cost of a decreased plant efficiency due to 
recirculation, D is the cost of damages imposed on the en- 
vironment, and C IS the cost of construction of the marine 
installations of the cooling water system. 

The term D is controversial and quantification in terms of 
money is elusive, if at all possible. Nevertheless this 
aspect of the waste heat problem plays a major role in plant 
site selection considerations- For power plants located at 
rivers this is understandable since waste heat disposal can 
cause temperature barriers, whereas it is less obvious for 
coastal locations where the influence of the waste heat 
disposal IS of a local character. 

Another characteristic feature of the term D is that there 
are strict limits as to the extent by which the environmen- 
tal impact can be managed by the lay-out of the cooling 
water system for a given plant site and a given rate of 
heat rejection. 

The term R which represents the power loss due to recircula- 
tion can be influenced to a great extent by a proper deter- 
mination of the positioning of intake and outfall. 


Recirculation Optimization 

It appears worthwhile to consider the objective function 


F = R + C (2) 

to investigate the consequences of this approach, and to 
raise the question of a possible conflict between a recir- 
culation optimization and a total optimization (equation 
( 1 ) ) . 

This approach is clearly applicable to cooling ponds where 
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environmental damages to the equatic environment need not 
be considered. With the numerical modeling techniques 
available today the elements necessary to perform a recir- 
culation optimization are available. 

Differentiation of the reduced objective function leads to 
the simple result that a total cost minimum is found when 
the marginal costs of decreasing the intake temperature 
equals the marginal cost of efficiency loss due to recircu- 
lation; 


dR _ dC 

dATi dATi (3) 

where ATi is the water temperature above its natural equili- 
brium at the intake. 

It should be noted that apart from the hydraulic complexity 
of the problem, a rigorous analysis is difficult to perform 
since the cost of efficiency loss can be evaluated in 
numerous ways. Furthermore, the effect of a certain degree 
of recirculation can to some extent be catered for in the 
design of the steam condensers. Regardless of the difficul- 
ties inherent in the problem it seems appropriate to perform 
an analysis of this type rather than to rely on arbitrary 
criteria with respect to the acceptable level of recircula- 
tion . 

Note that the term recirculation designates indirect recir- 
culation (intake located in the far-field) . Direct passage 
of heated water from outfall to intake is assumed to be non- 
existent . 


PLANT SITE SELECTION 

In principle, all terms in the objective function should be 
evaluated to provide a proper basis for a decision with re- 
spect to selection of potential sites. This in turn would 
require a full investigation of hydraulic and biological 
aspects at each site. However, a general knowledge with re- 
spect to flow patterns of the site is often sufficient as 
a basis for preliminary excess temperature field computa- 
tions and subsequent evaluation of biological effects. In 
this way the various sites can be given priority, and a few 
be selected for detailed hydraulic and biological studies. 

In the following, the main features of a site selection are 
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discussed, and a sample computatxon of the optimum recircu- 
lation IS given. 


Total Energy Budget 

The ultimate sink of waste heat is the atmosphere. Under 
equilibrium conditions, or as an average over a sufficient- 
ly long period of time, the total rate of heat dissipation 
over the influenced water surface area therefore must equal 
the rate of heat supply. 

For simplicity the rate of additional heat loss which takes 
place to the atmosphere when the temperature is raised AT 
above its natural equilibrium is assumed to be given by; 

H = f AT Watt/m2 

Conservation of energy requires that the surface rate of 
heat loss integrated over the entire excess temperature 
field equals the rate of heat supply; 

QATop Cp = f/ATdA (5) 

A 

where Q is the cooling water discharge, ATo is the tempe- 
rature increase over the condensers, p is the density of 
water and Cp is the specific heat. 

Evidently, the magnitude of the integral in "equation (5) is 
fixed for a given heat supply and completely independent of 
the currents , the mixing processes of the receiving water 
body and the mode of discharge. Furthermore, the integral 
IS independent of the vertical distribution of the excess 
temperature. These factors only affect the shape of the tem- 
perature field. Thus, in discharge areas with a considerable 
tidal activity the excess temperature distribution will tend 
to attain a more uniform distribution that in areas where 
the natural energy available for the spreading is lower. 

The result will be that the excess temperatures in receiving 
waters with a high level of mixing in the vicinity of the 
outfall will be relatively low, whereas the excess tempera- 
tures further away will be relatively high. Thus, the area 
encircled by lines through equal low excess temperatures 
will increase with increasing tidal and mixing intensity of 
the receiving water. 

It IS worthwhile to note that the shape of the excess tem- 
perature field IS determined to a very high extent by the 
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ability of the tidal currents to yield an initial spreading 
of the source. This appears to be a factor of considerable 
importance as compared to the effect v/hich can be obtained 
by a submerged diffuser. 

If a high velocity surface jet is utilized the jet momentum 
can induce a current pattern which can be of considerable 
importance for the shape of the entire temperature field 
in situations where the tidal currents are weak. 

Contrary to the distribution of the surface excess tempera- 
tures, the total amount of heat stored in the water depends 
on the transport and mixing characteristics of the receiving 
water as well as the mode of discharge. The heat storage is 
reduced by floating the heated water onto the receiving 
water. In this way the surface heat loss may not be signi- 
ficantly altered but the heat storage can be reduced consi- 
derably and temperature changes at the bottom of the re- 
ceiving water avoided. 


The Relevance of Initial Dilution 


A considerable attention has previously been given to the 
question of the so-called near field, i.e. the zone in the 
vicinity of the outfall where outfall conditions signifi- 
cantly affect the temperature field. This is especially re- 
levant for small discharges as far as evaluation of the 
risk of direct recirculation is concerned. Hov/ever, for 
cooling water discharges from large power plants heat stored 
in the discharge area can easily attain a significant level 
whereby the entraining water is already heated. 

The statement above especially applies to a situation in 
which the heated water is discharged through a diffuser in 
a receiving water characterized by an oscillating tidal 
flow. In this case discharge may take place over a diffuser 
length of a few hundred meters while the ridal current pro- 
vides a spreading of the source equal to the tidal excursi- 
on which may be several kilometers. Furthermore, since the 
diffuser jets entrain water which is already heated the di- 
lution obtained could be illusory. 

If, on the other hand, the heated water is discharged as a 
norizontal surface jet with an initial velocity comparable 
to the tidal velocity the momentum of the jet will induce 
a circulation pattern which may have a considerable effect 
on the temperature distribution. The temperature field is 
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then stretched in the direction of the jet and thus distor- 
ted as compared to the diffuser situation. However, regard- 
less of the mode of discharge the surface heat loss integra- 
ted over the whole warm water surface remains the same . 

Thermal regulations which prescribe maximum surface tempera- 
ture differentials can on the grounds above be said to have 
little meaning for large cooling water discharges. In this 
connection it should be emphasized that this aspect is onl^ 
to a limited extent manageable by the design of the intake 
and outfall structures. 

In order to make reliable predictions of the temperature 
field in the vicinity of the outfall the practising engineer 
IS faced with the problem of predicting the effect of the 
interaction between the near- and the far-field. A computa- 
tion based on the assumption that water unaffected by the 
discharge is available for dilution will inevitably lead to 
erroneous results. 


The Effect of Tidal Currents and Mass Injection 

An idealized example in which heat is injected to a semi- 
infinite region of open and shallow v;ater is considered. 

It IS furthermore assumed that the excess temperature is 
uniformly distributed over the depth d. The effect of a ti- 
dal current becomes apparent when the water is considered 
to remain at rest and the heat source to move sinusoidally 
over the tidal excursion, J.M, Williams [1] . The heat can in 
this situation be assumed to take place over a line source 
with a distribution and a resulting excess temperature 
field shown in principle in Fig. 1. 

An approximate solution with respect to the temperature dis- 
tribution in the center line of the oscillating temperature 
field can be obtained by considering the one-dimensional 
conservation of heat equation: 

^ _ _Q_ ^ (6) 

dx^ edD dx dDpCp 

where the symbols not previously defined are: 

X IS the coordinate in the direction perpendicular to the 
coast, D IS the dispersion coefficient in the direction per- 
pendicular to the coast and e is an equivalent tidal excur- 
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Sion • 


Equation (6) ignores the effect of heat dispersion in the 
direction of the current, and will for tnis reason give 
conservative estimates of the excess temperatures. This will 
particularly be the case for the lowest excess temperatures. 
However, the longitudinal dispersion effect can to some 
extent be catered] for by using a tidal excursion which is 
somewhat greater than the true excursion. 

The soulution to equation (6) is: 

AT = ATm exp (ax) (7) 


where 


a 


1 _Q_ 

2 edD 



feMP 

pCpQ^ 


( 8 ) 


and ATm is the maximum temperature along the coast subject 
to the condition that the total rate of heat dissipation is 
equal to the rate of supply: 


ATm = ATo (9) 

ef 

where ATo is the excess temperature at the point of discharge 
and Q is the cooling water flow. 


For large cooling water discharges Q and a small tidal ex- 
cursion equation (8) can be approximated by: 


a 


f e 

pC Q 
P 


( 10 ) 


in which case equation (7) is modified to: 


AT = ATo exp 


(- 


fe 

pCpQ 


X) 


( 11 ) 


For conditions normally encountered, the approximation used 
to derive equation (11) becomes good when the discharge 
exceeds 50-100 mVs corresponding to 1000-2000 MW nuclear 
power plants. Thus, for plants of this size the injection 
of mass becomes the dominating heat transport mechanism in 
the oscillating flow situation considered above. Note that 
equation (11) represents the solution to equation (6) when 
omitting the diffusive term. 
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Clearly, the idealized example above contains too many 
approximations to provide a proper basis for a determination 
of the positioning of intake and outfall. On the other hand, 
a useful insight with respect to orders of magnitude and 
parameter sensitivity can be gained with a very limited 
effort . 

Another idealized example which is useful to consider re- 
fers to the situation in which heat is injected over a ver- 
tical line source in an infinite region of shallow water 
with horizontal and isotropic diffusion. The conservation 
of heat equation is in this case: 

+ (l 1 ^ - — f — iT = 0 (12) 

V 2rdDy r dr pC^dD 

where r is the radial coordinate with origin at the source. 

The solution to this general equation contains the modified 
Bessel function of second kind and of order 1 Q 

4 TidD 

For tne special case wnen Q equals 2irdD the solution be- 
comes : 


4T = 4T0 exp r ) ' (13) 

For very large flows the solution is approximated by: 

AT = ATo exp (- r^) (14) 

Wore that the temperature at the source in this case, where 
no tidal activity is assumed to be present, becomes equal 
to me temperature of the source itself for relatively low 
discnarges. For conditions normally encountered this is the 
case for discharges in the order of 20 m^/s and above. 


The Effect of a Residual Drift 


The idealized examples above will rarely apply without 
modifications to real situations. For a more refined ana- 


8 


V-A-37 


lysis, two-dimensional time-dependent numerical models ba- 
sed on flow patterns generated by hydrodynamic models should 
be applied. Techniques applied in numerical modeling of this 
type is described in a paper presented by G.S. Rodenhuis at 
this Conference [2] . 

One of the reasons for the need for numerical modeling is 
the inability of analytical methods to incorporate the 
effect of irregular flow patterns, irregular topographies 
and the combined effect of the mass injection effect and 
the important effect of a residual current. 

Clearly, the presence of a residual drift current is of 
decisive importance for the shape of the temperature field, 
and hence for biological assessments and positioning of in- 
take and outfall. A residual drift current of, say, 3 cm/s 
would be about one magnitude greater than the advection 
velocity caused by the mass injection considered in the os- 
cillating flow example above. This will clearly have the 
effect of transforming the temperature field into a rela- 
tively narrow plume along the coast in the direction of the 
drift current. Furthermore, since unaffected water is con- 
tinuously fed into the plume, a considerable decrease of 
maximum excess temperatures and a corresponding increase of 
temperatures further downstream will take place. 

The recognition of the effect of a drift current should be 
duely reflected in the planning of field measurements and 
the subsequent data processing. 

In areas where the wind is a dominating factor in the gene- 
ration of the flow pattern which therefore is of a stochas- 
tic nature, a deterministic model approach is less suitable. 
As a compromise, deterministic numerical models can be run 
over a period long enough to contain the characteristic 
statistical features of the natural conditions in order to 
produce results which can be interpreted accordingly. 

Another possibility is to perform a number of model runs 
using parameter-sets which combine into a known frequency 
of occurrence. 

Finally, pure stochastic models can be employed. This 
approach is rarely used in this context, even though it 
would appear probable that this would be appropriate both 
for purposes of recirculation management and for assessments 
of biological effects. The advantages with respect to eluci- 
dation of both aspects appear to be self-evident, although 
it is unconceivable to apply a stochastic model without re- 
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sort to some type of deterministic approach. 


RECIRCULATION MANAGEMENT - AN EXAMPLE. 

According to equation (7) the optimum recirculation is found 
by equating the marginal costs of power loss and construc- 
tion costs; 


dR ^ dC 
dAT^ dAT^ 


Cost of Power Loss Due to Recirculation 

Harlemann [3] indicates the following relation for the 
change in basic plant efficiency for a typical boiling water 
or pressurized water reactor plant: 

^ ( 16 ) 

E - T^-T^ 

where Tjj-T^ is the deviation in (*^F) from the design steam 
condensing'*' temperature, T 2 is the maximum operating tempera- 
ture of the plant, and T^^'^is the design steam condensing 
temperature. Equation (16) , which qualitatively is correct 
[3] , shows that there is a nearly linear relationship be- 
tween power loss and deviations from the design temperature. 
Consequently, the present value of power loss due to the 
presence of temperatures increased AT^^ over the natural 
equilibrium at the intake can be written as; 

R = AE f AT^ PVF (17) 

where AE is the power loss per degree increase in intake 
water temperature, f is the cost of IMIV year, PVF is the 
present value factor. 

These factors are tentatively evaluated as follows: 

AE = 0.0035E MW ®c“^ 

f = 20 $ (MWh)“l = 175,000 $ (MW year) 

PVF = 14 for a 5% discount rate and a pe- 
riod of 25 years. 
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For a 2000 MW power plant the marginal present value of power 
loss IS estimated at 17 million dollars per degree centigrade. 

Construction Costs . 

Jn the example considered above it is assumed that the 
plant's required 100 m^/s of cooling water is conveyed 
through large submerged conduits from an intake structure 
several hundred meters from the shore, see Fig. 1. The out- 
fall is assumed to be located at the shore. For simplicity 
It is furthermore assumed that the construction and pumping 
costs can be expressed as: 


C = Px (18) 

P is the cost per unit lenght and x is the distance between 
intake and outfall. 


The time averaged excess temperature distribution in the 
outfall-intake line is assumed approximated by an exponen- 
tial decay, equation (7) . The marginal cost of construction 
and ptamping can then be computed as: 


dC = _P 

dAT^ aAT^ 


(19) 


Substitution of equation (17) and (18) in (3) leads to the 
optimal excess temperature at the intake: 


AT optim = 

1 ^ AE f PVF a 

from which the optimal distance between outfall and intake 
is readily computed using equation (7) . 

The cost of the marine installations (including pumping) of 
the type considered is estimated at 80,000 US $ per meter. 

The exponent a in equation (7) is assumed to be increased 
by the action of a residual tidal drift current as compared 
to the value one would obtain in the case of a purely os- 
cillating flow. The magnitude of the exponent is estimated 
at 0.003 m“l. 

With the above cited figures the optimum recirculation be- 
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comes 


AT^optim 


80 10^ 

17 10® 3 10 


1.6 


o 


C. 


Because of the exponential decay of the excess 
temperature with distance from the outfall and because of 
the considerable construction costs involved, deviations 
from the optimum temperature will inevitably lead to a con- 
siderable increase in the total cost. 

The choice of magnitude of the various parameters involved 
in the analysis above can of course be discussed and a sen- 
sitivity test would in an actual case be required. Further- 
more, the effect of increased ambient temperatures in the 
far-field can to some extent be catered for in the design 
of the steaim condensers. However, in any case, the analysis 
stresses the importance of a reliable prediction of the 
excess temperature distribution. 

As a practical comment to the example considered it should 
be noted that it would appear to be advantageous to dis- 
charge in the off-shore position using a high velocity jet, 
and to take water in close to the shore. Furthermore, a 
position of the intake upstream relative to the outfall con 
duits should be considered. Thus, in order to perform a 
complete analysis, two dimensions for determination of the 
optimum intake position would need be considered. 


RECIRCULATION AND BIOLOGICAL EFFECTS 

The question whether or not decisions based on recircula- 
tion management principles could be in conflict with biolo- 
gical demands is difficult to assess. Apart from the very 
limited effect of the efficiency decrease due to recircula- 
tion, the same quantity of heat is rejected to the environ- 
ment. The fact that recirculation is present only means 
that the excess temperature field is "concentrated" to a 
higher extent in the area between outfall and intake than 
would be the case if no recirculation was present at all. 

Additionally, it appears that a recirculation management 
analysis would always tend to yield solutions which are pre 
ferable from a biological point of view as well since sites 
and intake-outfall lay-out solutions causing relatively 
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small areas to be subject to high excess temperatures would 
be preferred for both reasons. 

There may , however , be coastal sites which for biological 
reasons should not be used. It should be noted that as far 
as decimation of marine organism due to entrainment in the 
cooling water system is concerned, this aspect is indepen- 
dent of the question of the distribution of an excess tem- 
perature field. This and other biological aspects can be 
assessed using a numerical modeling technique described in 
the paper presented by B. Miller and K.I. Dahl-Madsen lat 
this Conference [4]. 


DISCUSSION 

The aim of the considerations presented in this paper has 
been to elucidate some of the important elements in the 
analysis of the was.te heat management problem. Therefore, 
the methods indicated above do by no means represenr a 
state-of-the-art of modeling techniques but should be regar- 
ded as tools which can be used for purposes of quick appro- 
ximations, sensitivity analyses and systems analysis. 

The main object of this paper is that it addresses the 
genuine need to recognize that the waste heat problem is 
not simply solved by prescribing maximum surface tempera- 
ture differentials in the receiving water and requiring zero 
recirculation. 
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IDEALIZED TEMPERATURE DISTRIBUTIONS AT THE SLACK 
TIDES OF AN OSCILLATING TIDAL FLOW FIELD 



ILLUSTRATIONS OF THE RECIRCULATION MANAGEMENT EXAMPLE 

laT, *c, excess reMPeRATuRe 


aT, oot r 


xopt 


METERS FROM OUTFALL 



AT, opt 5 aT, *C excess temperature 


Fig 1 Idealized temperature distribution and marginal 
cost analysis 


